Expression characteristics and potential function of non-coding RNA in mouse cortical cells.

Non-coding RNAs (ncRNAs) play essential regulatory functions in various physiological and pathological processes in the brain. To systematically characterize the ncRNA profile in cortical cells, we downloaded single-cell SMART-Seq v4 data of mouse cerebral cortex. Our results revealed that the ncRNAs alone are sufficient to define the identity of most cortical cell types. We identified 1,600 ncRNAs that exhibited cell type specificity, even yielding to distinguish microglia from perivascular macrophages with ncRNA. Moreover, we characterized cortical layer and region specific ncRNAs, in line with the results by spatial transcriptome (ST) data. By constructing a co-expression network of ncRNAs and protein-coding genes, we predicted the function of ncRNAs. By integrating with genome-wide association studies data, we established associations between cell type-specific ncRNAs and traits related to neurological disorders. Collectively, our study identified differentially expressed ncRNAs at multiple levels and provided the valuable resource to explore the functions and dysfunctions of ncRNAs in cortical cells.

Multimodal cell atlas of the ageing human skeletal muscle.

Muscle atrophy and functional decline (sarcopenia) are common manifestations of frailty and are critical contributors to morbidity and mortality in older people1. Deciphering the molecular mechanisms underlying sarcopenia has major implications for understanding human ageing2. Yet, progress has been slow, partly due to the difficulties of characterizing skeletal muscle niche heterogeneity (whereby myofibres are the most abundant) and obtaining well-characterized human samples3,4. Here we generate a single-cell/single-nucleus transcriptomic and chromatin accessibility map of human limb skeletal muscles encompassing over 387,000 cells/nuclei from individuals aged 15 to 99 years with distinct fitness and frailty levels. We describe how cell populations change during ageing, including the emergence of new populations in older people, and the cell-specific and multicellular network features (at the transcriptomic and epigenetic levels) associated with these changes. On the basis of cross-comparison with genetic data, we also identify key elements of chromatin architecture that mark susceptibility to sarcopenia. Our study provides a basis for identifying targets in the skeletal muscle that are amenable to medical, pharmacological and lifestyle interventions in late life.

A spatiotemporal atlas of cholestatic injury and repair in mice.

Cholestatic liver injuries, characterized by regional damage around the bile ductular region, lack curative therapies and cause considerable mortality. Here we generated a high-definition spatiotemporal atlas of gene expression during cholestatic injury and repair in mice by integrating spatial enhanced resolution omics sequencing and single-cell transcriptomics. Spatiotemporal analyses revealed a key role of cholangiocyte-driven signaling correlating with the periportal damage-repair response. Cholangiocytes express genes related to recruitment and differentiation of lipid-associated macrophages, which generate feedback signals enhancing ductular reaction. Moreover, cholangiocytes express high TGFß in association with the conversion of liver progenitor-like cells into cholangiocytes during injury and the dampened proliferation of periportal hepatocytes during recovery. Notably, Atoh8 restricts hepatocyte proliferation during 3,5-diethoxycarbonyl-1,4-dihydro-collidin damage and is quickly downregulated after injury withdrawal, allowing hepatocytes to respond to growth signals. Our findings lay a keystone for in-depth studies of cellular dynamics and molecular mechanisms of cholestatic injuries, which may further develop into therapies for cholangiopathies.

Region-specific transcriptomic responses to obesity and diabetes in macaque hypothalamus.

The hypothalamus plays a crucial role in the progression of obesity and diabetes; however, its structural complexity and cellular heterogeneity impede targeted treatments. Here, we profiled the single-cell and spatial transcriptome of the hypothalamus in obese and sporadic type 2 diabetic macaques, revealing primate-specific distributions of clusters and genes as well as spatial region, cell-type-, and gene-feature-specific changes. The infundibular (INF) and paraventricular nuclei (PVN) are most susceptible to metabolic disruption, with the PVN being more sensitive to diabetes. In the INF, obesity results in reduced synaptic plasticity and energy sensing capability, whereas diabetes involves molecular reprogramming associated with impaired tanycytic barriers, activated microglia, and neuronal inflammatory response. In the PVN, cellular metabolism and neural activity are suppressed in diabetic macaques. Spatial transcriptomic data reveal microglia's preference for the parenchyma over the third ventricle in diabetes. Our findings provide a comprehensive view of molecular changes associated with obesity and diabetes.

Uncovering the Bronchoalveolar Single-Cell Landscape of Patients With Pulmonary Tuberculosis With Human Immunodeficiency Virus Type 1 Coinfection.

BACKGROUND: Coinfection of human immunodeficiency virus type 1 (HIV-1) is the most significant risk factor for tuberculosis (TB). The immune responses of the lung are essential to restrict the growth of Mycobacterium tuberculosis and avoid the emergence of the disease. Nevertheless, there is still limited knowledge about the local immune response in people with HIV-1-TB coinfection. METHODS: We employed single-cell RNA sequencing (scRNA-seq) on bronchoalveolar lavage fluid from 9 individuals with HIV-1-TB coinfection and 10 with pulmonary TB. RESULTS: A total of 19 058 cells were grouped into 4 major cell types: myeloid cells, T/natural killer (NK) cells, B cells, and epithelial cells. The myeloid cells and T/NK cells were further divided into 10 and 11 subsets, respectively. The proportions of dendritic cell subsets, CD4+ T cells, and NK cells were lower in the HIV-1-TB coinfection group compared to the TB group, while the frequency of CD8+ T cells was higher. Additionally, we identified numerous differentially expressed genes between the CD4+ and CD8+ T-cell subsets between the 2 groups. CONCLUSIONS: HIV-1 infection not only affects the abundance of immune cells in the lungs but also alters their functions in patients with pulmonary TB.

VGLL1 cooperates with TEAD4 to control human trophectoderm lineage specification.

In contrast to rodents, the mechanisms underlying human trophectoderm and early placenta specification are understudied due to ethical barriers and the scarcity of embryos. Recent reports have shown that human pluripotent stem cells (PSCs) can differentiate into trophectoderm (TE)-like cells (TELCs) and trophoblast stem cells (TSCs), offering a valuable in vitro model to study early placenta specification. Here, we demonstrate that the VGLL1 (vestigial-like family member 1), which is highly expressed during human and non-human primate TE specification in vivo but is negligibly expressed in mouse, is a critical regulator of cell fate determination and self-renewal in human TELCs and TSCs derived from naïve PSCs. Mechanistically, VGLL1 partners with the transcription factor TEAD4 (TEA domain transcription factor 4) to regulate chromatin accessibility at target gene loci through histone acetylation and acts in cooperation with GATA3 and TFAP2C. Our work is relevant to understand primate early embryogenesis and how it differs from other mammalian species.

Presence of onco-fetal neighborhoods in hepatocellular carcinoma is associated with relapse and response to immunotherapy.

Onco-fetal reprogramming of the tumor ecosystem induces fetal developmental signatures in the tumor microenvironment, leading to immunosuppressive features. Here, we employed single-cell RNA sequencing, spatial transcriptomics and bulk RNA sequencing to delineate specific cell subsets involved in hepatocellular carcinoma (HCC) relapse and response to immunotherapy. We identified POSTN+ extracellular matrix cancer-associated fibroblasts (EM CAFs) as a prominent onco-fetal interacting hub, promoting tumor progression. Cell-cell communication and spatial transcriptomics analysis revealed crosstalk and co-localization of onco-fetal cells, including POSTN+ CAFs, FOLR2+ macrophages and PLVAP+ endothelial cells. Further analyses suggest an association between onco-fetal reprogramming and epithelial-mesenchymal transition (EMT), tumor cell proliferation and recruitment of Treg cells, ultimately influencing early relapse and response to immunotherapy. In summary, our study identifies POSTN+ CAFs as part of the HCC onco-fetal niche and highlights its potential influence in EMT, relapse and immunotherapy response, paving the way for the use of onco-fetal signatures for therapeutic stratification.

Unsupervised spatially embedded deep representation of spatial transcriptomics.

Optimal integration of transcriptomics data and associated spatial information is essential towards fully exploiting spatial transcriptomics to dissect tissue heterogeneity and map out inter-cellular communications. We present SEDR, which uses a deep autoencoder coupled with a masked self-supervised learning mechanism to construct a low-dimensional latent representation of gene expression, which is then simultaneously embedded with the corresponding spatial information through a variational graph autoencoder. SEDR achieved higher clustering performance on manually annotated 10 × Visium datasets and better scalability on high-resolution spatial transcriptomics datasets than existing methods. Additionally, we show SEDR's ability to impute and denoise gene expression (URL: https://github.com/JinmiaoChenLab/SEDR/ ).

STOmicsDB: a comprehensive database for spatial transcriptomics data sharing, analysis and visualization.

Recent technological developments in spatial transcriptomics allow researchers to measure gene expression of cells and their spatial locations at the single-cell level, generating detailed biological insight into biological processes. A comprehensive database could facilitate the sharing of spatial transcriptomic data and streamline the data acquisition process for researchers. Here, we present the Spatial TranscriptOmics DataBase (STOmicsDB), a database that serves as a one-stop hub for spatial transcriptomics. STOmicsDB integrates 218 manually curated datasets representing 17 species. We annotated cell types, identified spatial regions and genes, and performed cell-cell interaction analysis for these datasets. STOmicsDB features a user-friendly interface for the rapid visualization of millions of cells. To further facilitate the reusability and interoperability of spatial transcriptomic data, we developed standards for spatial transcriptomic data archiving and constructed a spatial transcriptomic data archiving system. Additionally, we offer a distinctive capability of customizing dedicated sub-databases in STOmicsDB for researchers, assisting them in visualizing their spatial transcriptomic analyses. We believe that STOmicsDB could contribute to research insights in the spatial transcriptomics field, including data archiving, sharing, visualization and analysis. STOmicsDB is freely accessible at https://db.cngb.org/stomics/.

Single cell multi-omics reveal intra-cell-line heterogeneity across human cancer cell lines.

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

Live birth of chimeric monkey with high contribution from embryonic stem cells.

A formal demonstration that mammalian pluripotent stem cells possess preimplantation embryonic cell-like (naive) pluripotency is the generation of chimeric animals through early embryo complementation with homologous cells. Whereas such naive pluripotency has been well demonstrated in rodents, poor chimerism has been achieved in other species including non-human primates due to the inability of the donor cells to match the developmental state of the host embryos. Here, we have systematically tested various culture conditions for establishing monkey naive embryonic stem cells and optimized the procedures for chimeric embryo culture. This approach generated an aborted fetus and a live chimeric monkey with high donor cell contribution. A stringent characterization pipeline demonstrated that donor cells efficiently (up to 90%) incorporated into various tissues (including the gonads and placenta) of the chimeric monkeys. Our results have major implications for the study of primate naive pluripotency and genetic engineering of non-human primates.

GDF11 slows excitatory neuronal senescence and brain ageing by repressing p21.

As a major neuron type in the brain, the excitatory neuron (EN) regulates the lifespan in C. elegans. How the EN acquires senescence, however, is unknown. Here, we show that growth differentiation factor 11 (GDF11) is predominantly expressed in the EN in the adult mouse, marmoset and human brain. In mice, selective knock-out of GDF11 in the post-mitotic EN shapes the brain ageing-related transcriptional profile, induces EN senescence and hyperexcitability, prunes their dendrites, impedes their synaptic input, impairs object recognition memory and shortens the lifespan, establishing a functional link between GDF11, brain ageing and cognition. In vitro GDF11 deletion causes cellular senescence in Neuro-2a cells. Mechanistically, GDF11 deletion induces neuronal senescence via Smad2-induced transcription of the pro-senescence factor p21. This work indicates that endogenous GDF11 acts as a brake on EN senescence and brain ageing.

Creation of a High-Throughput Microfluidic Platform for Single-Cell Transcriptome Sequencing of Cell-Cell Interactions.

Cell-cell interaction is one of the major modalities for transmitting information between cells and activating the effects of functional cells. However, the construction of high-throughput analysis technologies from cell omics focusing on the impact of interactions of functional cells on targets has been relatively unexplored. Here, they propose a droplet-based microfluidic platform for cell-cell interaction sequencing (c-c-seq) and screening in vitro to address this challenge. A class of interacting cells is pre-labeled using cell molecular tags, and additional single-cell sequencing reagents are introduced to quickly form functional droplet mixes. Lastly, gene expression analysis is used to deduce the impact of the interaction, while molecular sequence tracing identifies the type of interaction. Research into the active effect between antigen-presenting cells and T cells, one of the most common cell-to-cell interactions, is crucial for the advancement of cancer therapy, particularly T cell receptor-engineered T cell therapy. As it allows for high throughput, this platform is superior to well plates as a research platform for cell-to-cell interactions. When combined with the next generation of sequencing, the platform may be able to more accurately evaluate interactions between epitopes and receptors and verify their functional relevance.

A scATAC-seq atlas of chromatin accessibility in axolotl brain regions.

Axolotl (Ambystoma mexicanum) is an excellent model for investigating regeneration, the interaction between regenerative and developmental processes, comparative genomics, and evolution. The brain, which serves as the material basis of consciousness, learning, memory, and behavior, is the most complex and advanced organ in axolotl. The modulation of transcription factors is a crucial aspect in determining the function of diverse regions within the brain. There is, however, no comprehensive understanding of the gene regulatory network of axolotl brain regions. Here, we utilized single-cell ATAC sequencing to generate the chromatin accessibility landscapes of 81,199 cells from the olfactory bulb, telencephalon, diencephalon and mesencephalon, hypothalamus and pituitary, and the rhombencephalon. Based on these data, we identified key transcription factors specific to distinct cell types and compared cell type functions across brain regions. Our results provide a foundation for comprehensive analysis of gene regulatory programs, which are valuable for future studies of axolotl brain development, regeneration, and evolution, as well as on the mechanisms underlying cell-type diversity in vertebrate brains.

Single-cell chromatin accessibility profiling of cell-state-specific gene regulatory programs during mouse organogenesis.

In mammals, early organogenesis begins soon after gastrulation, accompanied by specification of various type of progenitor/precusor cells. In order to reveal dynamic chromatin landscape of precursor cells and decipher the underlying molecular mechanism driving early mouse organogenesis, we performed single-cell ATAC-seq of E8.5-E10.5 mouse embryos. We profiled a total of 101,599 single cells and identified 41 specific cell types at these stages. Besides, by performing integrated analysis of scATAC-seq and public scRNA-seq data, we identified the critical cis-regulatory elements and key transcription factors which drving development of spinal cord and somitogenesis. Furthermore, we intersected accessible peaks with human diseases/traits-related loci and found potential clinical associated single nucleotide variants (SNPs). Overall, our work provides a fundamental source for understanding cell fate determination and revealing the underlying mechanism during postimplantation embryonic development, and expand our knowledge of pathology for human developmental malformations.

Single-nucleus chromatin landscapes during zebrafish early embryogenesis.

Vertebrate embryogenesis is a remarkable process, during which numerous cell types of different lineages arise within a short time frame. An overwhelming challenge to understand this process is the lack of dynamic chromatin accessibility information to correlate cis-regulatory elements (CREs) and gene expression within the hierarchy of cell fate decisions. Here, we employed single-nucleus ATAC-seq to generate a chromatin accessibility dataset on the first day of zebrafish embryogenesis, including 3.3 hpf, 5.25 hpf, 6 hpf, 10 hpf, 12 hpf, 18 hpf and 24 hpf, obtained 51,620 high-quality nuclei and 23 clusters. Furthermore, by integrating snATAC-seq data with single-cell RNA-seq data, we described the dynamics of chromatin accessibility and gene expression across developmental time points, which validates the accuracy of the chromatin landscape data. Together, our data could serve as a fundamental resource for revealing the epigenetic regulatory mechanisms of zebrafish embryogenesis.

High-Throughput Screening of Functional Neo-Antigens and Their Specific T-Cell Receptors via the Jurkat Reporter System Combined with Droplet Microfluidics.

T-cell receptor (TCR)-engineered T cells can precisely recognize a broad repertoire of targets derived from both intracellular and surface proteins of tumor cells. TCR-T adoptive cell therapy has shown safety and promising efficacy in solid tumor immunotherapy. However, antigen-specific functional TCR screening is time-consuming and expensive, which limits its application clinically. Here, we developed a novel integrated antigen-TCR screening platform based on droplet microfluidic technology, enabling high-throughput peptide-major histocompatibility complex (pMHC)-to-TCR paired screening with a high sensitivity and low background signal. We introduced DNA barcoding technology to label peptide antigen candidate-loaded antigen-presenting cells and Jurkat reporter cells to check the specificity of pMHC-TCR candidates. Coupled with the next-generation sequencing pipeline, interpretation of the DNA barcodes and the gene expression level of the Jurkat T-cell activation pathway provided a clear peptide-MHC-TCR recognition relationship. Our proof-of-principle study demonstrates that the platform could achieve pMHC-TCR paired high-throughput screening, which is expected to be used in the cross-reactivity and off-target high-throughput paired testing of candidate pMHC-TCRs in clinical applications.

An invasive zone in human liver cancer identified by Stereo-seq promotes hepatocyte-tumor cell crosstalk, local immunosuppression and tumor progression.

Dissecting and understanding the cancer ecosystem, especially that around the tumor margins, which have strong implications for tumor cell infiltration and invasion, are essential for exploring the mechanisms of tumor metastasis and developing effective new treatments. Using a novel tumor border scanning and digitization model enabled by nanoscale resolution-SpaTial Enhanced REsolution Omics-sequencing (Stereo-seq), we identified a 500 µm-wide zone centered around the tumor border in patients with liver cancer, referred to as "the invasive zone". We detected strong immunosuppression, metabolic reprogramming, and severely damaged hepatocytes in this zone. We also identified a subpopulation of damaged hepatocytes with increased expression of serum amyloid A1 and A2 (referred to collectively as SAAs) located close to the border on the paratumor side. Overexpression of CXCL6 in adjacent malignant cells could induce activation of the JAK-STAT3 pathway in nearby hepatocytes, which subsequently caused SAAs' overexpression in these hepatocytes. Furthermore, overexpression and secretion of SAAs by hepatocytes in the invasive zone could lead to the recruitment of macrophages and M2 polarization, further promoting local immunosuppression, potentially resulting in tumor progression. Clinical association analysis in additional five independent cohorts of patients with primary and secondary liver cancer (n = 423) showed that patients with overexpression of SAAs in the invasive zone had a worse prognosis. Further in vivo experiments using mouse liver tumor models in situ confirmed that the knockdown of genes encoding SAAs in hepatocytes decreased macrophage accumulation around the tumor border and delayed tumor growth. The identification and characterization of a novel invasive zone in human cancer patients not only add an important layer of understanding regarding the mechanisms of tumor invasion and metastasis, but may also pave the way for developing novel therapeutic strategies for advanced liver cancer and other solid tumors.

Spatially informed clustering, integration, and deconvolution of spatial transcriptomics with GraphST.

Spatial transcriptomics technologies generate gene expression profiles with spatial context, requiring spatially informed analysis tools for three key tasks, spatial clustering, multisample integration, and cell-type deconvolution. We present GraphST, a graph self-supervised contrastive learning method that fully exploits spatial transcriptomics data to outperform existing methods. It combines graph neural networks with self-supervised contrastive learning to learn informative and discriminative spot representations by minimizing the embedding distance between spatially adjacent spots and vice versa. We demonstrated GraphST on multiple tissue types and technology platforms. GraphST achieved 10% higher clustering accuracy and better delineated fine-grained tissue structures in brain and embryo tissues. GraphST is also the only method that can jointly analyze multiple tissue slices in vertical or horizontal integration while correcting batch effects. Lastly, GraphST demonstrated superior cell-type deconvolution to capture spatial niches like lymph node germinal centers and exhausted tumor infiltrating T cells in breast tumor tissue.

Cell atlas of CCl 4-induced progressive liver fibrosis reveals stage-specific responses.

Chronic liver injury leads to progressive liver fibrosis and ultimately cirrhosis, a major cause of morbidity and mortality worldwide. However, there are currently no effective anti-fibrotic therapies available, especially for late-stage patients, which is partly attributed to the major knowledge gap regarding liver cell heterogeneity and cell-specific responses in different fibrosis stages. To reveal the multicellular networks regulating mammalian liver fibrosis from mild to severe phenotypes, we generated a single-nucleus transcriptomic atlas encompassing 49 919 nuclei corresponding to all main liver cell types at different stages of murine carbon tetrachloride (CCl 4)-induced progressive liver fibrosis. Integrative analysis distinguished the sequential responses to injury of hepatocytes, hepatic stellate cells and endothelial cells. Moreover, we reconstructed cell-cell interactions and gene regulatory networks implicated in these processes. These integrative analyses uncovered previously overlooked aspects of hepatocyte proliferation exhaustion and disrupted pericentral metabolic functions, dysfunction for clearance by apoptosis of activated hepatic stellate cells, accumulation of pro-fibrotic signals, and the switch from an anti-angiogenic to a pro-angiogenic program during CCl 4-induced progressive liver fibrosis. Our dataset thus constitutes a useful resource for understanding the molecular basis of progressive liver fibrosis using a relevant animal model.